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\ HST NICMOS narrowband images of the shocked molecular hydrogen emis- 

^ \ sion in OMC-1 are analyzed to reveal new information on the BN/KL outflow. 

I— 1 1 The outstanding morphological feature of this region is the array of molecular 

"■q^I hydrogen "fingers" emanating from the general vicinity of IRc2 and the presence 

Q \ of several Herbig-Haro objects. The NICMOS images appear to resolve individ- 

^ \ ual shock fronts. This work is a more quantitative and detailed analysis of our 

\ data from a previous paper (Schultz et al.). 

\ Line strengths for the H2 1-0 S(4) plus 2-1 S(6) lines at 1.89 /xm are estimated 

^ I from measurements with the Paschen a continuum filter F190N at 1.90 /um, and 

Q^ \ continuum measurements at 1.66 and 2.15 /im. We compare the observed H2 line 

I strengths and ratios of the 1.89 /im and 2.12 /im 1-0 S(l) lines with models for 

CN \ molecular cloud shock waves. Most of the data cannot be fit by J-shocks, but 

0\ ' are well matched by C-shocks with shock velocities in the range of 20-45 km s~^ 

' and preshock densities of 10^ — 10^ cm~'^, similiar to values obtained in larger 

^ ' beam studies which averaged over many shocks. There is also some evidence that 

^ ' shocks with higher densities have lower velocities. 

I Subject headings: H II regions - infrared: ISM: lines and bands - ISM: individual 

(OMC-1) - ISM: jets and outflows - stars:pre- main-sequence 
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Introduction 



At 450 pc, the Orion molecular cloud is the nearest and best-studied region of massive 
star formation. The Trapezium stars, formed within Orion Molecular Cloud 1, have cleared 
a cavity at the near edge of the cloud. The visible Orion Nebula is the thin layer of photo- 
ionized gas on the cavity's surface facing the observer. Behind M42, and further from the 
observer, is a photodissociation region (PDR) also excited by the Trapezium. The BN/KL 
region lies still deeper in the cloud, beyond both the ionized gas and the PDR. This region 
contains embedded sources wit h one or more associate d outflows; the total luminosity of 
this region approaches ~ IO^Lq (IGenzel fc Stutzkilll989l ). The mid- infrared s ource IRc2 had 
long been thought to be the origin of these outflows; but iDougados et al.l (Il993l ) resolved 
IRc2 into four sources, raisi ng the possibility that none of them is sufficiently powerful to 
drive the observed outflows . iMenten fc Reid I (I19951) suggested that the origi n may be closer 
to the infrared source "n" (ILonsdale et al.lll982l ). located ~ 5" SW of IRc2. iGreenhill et al. 



( 120041 ) detected extended emission from source n at wavelengths out to 22 /im, but estimate 
a luminosity of only 2000 Lq. They also resolved IRc2 into ~ 5 knots and suggested that 
these sources together with radio source I (IChurchwell et al.lll987l ) comprise at least part of 
the core of a high density star forming cluster. 

The most striking of the molecular hydrogen outflows is an ~ 3' (0.4 pc) sized, butterfly- 
shaped region of H2 emission, centered to the north of BN, w hich exhibits line ratio s typical of 
shock excitation (IBeckwith et al.lll978l ). From O I emission. I Axon fc Tavlorl (119841 ) identified 



a number of optical HH objects in this vicinity. 



H2 structures in the outfiow. lAUen fc BurtonI (Il993l . hereafter AB) showed that these H2 



Tavlor et al.l ( 119841 ) discovered peculiar linear 



"fingers" and all the associated optical HH objects at the far norther n end of the outfiow 
(approximately 120" from BN) terminated in knots of Fe II emission. IStolovy et al.l (119981 ) 
found additional H2 fingers within 3 0" of BN. Schultz et al. (1999, Paper 1), found that 
only 2 of the 15 inner fingers seen by IStolovy et al.l (119981 ) had bow shocks capped by Fe II 
emission, suggesting a lower excitation than in the outer fingers seen by AB. 



From offsets between the peak H2 emission and the peak H2 velocity, iGustafsson et al. 



(I2OO3I ) suggested that the H2 emission arises in part from outflows from protostars within 
dense clu r nps o f gas. In contrast, based on proper motion studies of optical features, 
Doi et al.l (120021 ) found that both flnger systems could have been created by an explosive 
event close to the IRc2-BN comp lex w hich took pla c e app roximately 1000 years ago. In- 
terestingly, iRodriguez et al.l (120051 ) and lGomez et al.l (120051 ) suggested that BN and sources 
I and n were originally part of a multiple massive stellar system that disintegrated about 
500 years ago. Explanations for the unique system of flngers have focused on two theories. 
AB originally suggested that they are "bullets" — ejected clumps leaving a wake of shocked 
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material behind them. Howe ver, the observed morphology of the H2 e r aission is inconsis- 



tent with inodels for bullets (IStone fc NormanI Il992l . iKlein et al.l 1 1994 IXu &: Stond Il995 



Jones et al.l Il996l ). These models predict that rapidly moving clumps are fragmented and 
also p redict that the tai ls should be pointing away from the ejection source, which is not 



seen. 



Stone et al.l (Il995l ) suggested the features are produced when a faster wind collides 



with a slower, older outflow. Rayleigh- Taylor instabilities from the collision form the clumps 
in situ, moving at the speed of the older outflow. The observed fingers then condense be- 
hind the slowly moving clumps. This is si milar to the mechanism t hought to have produced 
the cometary knots in the Helix Nebula (lO'Dell fc Handrorull996l ). One prediction of this 
model is that a region of cl u mpy H^ emission will form behind (i.e. upstream of) the bullets. 
McCaughrean fc Mac Lowl (119971 ) claimed to have foun d this clu mpy emission in the central 



regio n of the H2 outflow. However, our previous work (jPaper ll ) and that of (IStolovy et al. 



19981 ) shows that much of this "clumpy" H2 emission is resolved into more discrete objects, 
some resembling additional fingers. The remaining, unresolved clumpy emission is often 
mixed with the inner fingers. 

Here we discuss the interpretation of our previously published NICMOS infrared images 
of a 90" wide region centered on BN/KL, focusing on the structure of the H2 emission. 
Examination of the F190N im ages, originally obtained for subtracting continuum from the 
Pa 1.87 /um images (jPaper 1 ). suggested that in many regions there was a strong correlation 
with the H2 1-0 S(l) 2.12 //m continuum-subtracted images. In fact, the F190N filter 
bandpass includes both the H^ 1-0 S(4) arid the 2-1 S(6) lines at 1.89 /im. This H2 emission 
is likely produced in shocks (iGautier et al.lll976l ). H2 emission can also be produced throug h 
UV fluorescence , but l arger beam studies of multiple H2 transitions by lUsuda et al.l (119961 ) . 
Rosenthal et al.l (120001 ). and many others have shown that the measured line ratios in this 
region are consistent only with thermal excitation and not UV fluorescence. In this paper, 
we focus on the finger-like structures and the HH objects. Some of these objects have optical 
counterparts, which places constraints on their position within the cloud/nebula interface. 
In §2 we discuss the observations and data reduction. In §3 we compare the observed H2 line 
brightnesses with several classes of shock models in order to determine shock types, shock 
velocities, and gas densities. In §4 we discuss the morphology of, and emission from, many 
of the more distinct, brighter features seen in our images. 



2. Observations 



Observations were made of H2 1-0 S(l) 2.12 fim and 1-0 S(4) plus 2-1 S(6) 1.89 /im in 
the 1% bandpass NICMOS filters. The initial reduction of the NICMOS data was described 
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in detail in Paper 1 and generally followed standard proced ures. Photom etry was then per- 
formed on the reduced data from Paper 1 using the IRAF (jTodylll993l ) task polyphot, in 
which the average brightness is estimated inside a user-defined polygonal aperture. The aper- 
tures were designed to closely follow the outline of each object we identified. Sky subtraction 
was not performed with the polyphot task, but separately using a region or regions far from 
areas of obvious emission. Regions selected for photometry are shown in Figure 1. Knot 



designations, except for HH 208, utilize the source identification scheme of lO 'Dell fc Wen 



(119941 ). Knot identifications for HH 208 are shown in Figure 11. The features were selected 
for being distinct and fairly bright, having detectable 1.66 /im and 2.15 ^m continuum emis- 
sion, and 2.12 fim line emission. Knot U was also included even though it had no detectable 
1.66 /im continuum. The resulting photometry is listed in Columns 2-5 of Table 1. The 
formal statistical errors are such that the signal-to-noise ratios of all the measurements in 
Table 1 exceed 25 except for three: i) the 1.66 fim measurement of 128-248 (S/N = 4), ii) 
the 1.66 /zm measurement of HH 208U (not detected), and iii) the 2.15 fim measurement of 
HH 208U (S/N=15). 

There are no continuum observations for the H2 1.892 fim 1-0 S(4) plus 2-1 S(6) images. 
However, we did observe continua at 1.66 and 2.15 /im, intended as continua for the Fe II and 
1-0 S(l) lines respectively. The 1.89 //m continuum has been estimated from these contin- 
uum measurements by linearly interpolating between the 1.66 /tm and 2.15 /tm photometric 
measurements. The continuum interpolation approach was checked by applying the same 
technique to nine regions with no 2.12 fim H2 emission. In these test regions, the interpo- 
lated continuum value was on average 99 ± 3% of the value actually measured in the F190N 
filter, confirming that this is a reasonable approach. As an example of the interpolation, the 
brightnesses for 128-248, with the estimated continuum is shown in Figure 2. 
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Fig. 1. — H2 1-0 S(l) image of the Orion region with the selected features labelled. North 
is up and east is to the left. BN is located at 05^35^14.^2, -05°22'23."2 (J2000) and the 
location of IRc2 is noted by a cross. 
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Table 1. H2 Photometry and Results 



Measured Line plus Continuum Extinction Corrected H2 Brightnesses 

Position 1.66 /im 1.90 fim 2.12^m 2.15 ^2.12^111'" Ay^ 1-89 fiin 2.12 

[10^^ ergs s^^ cm^^ sr^^] (mag) (niag) [10^^ ergs s^^ cm^^ sr^"*^] 



HH 208A 


0.31 


1.14 


2.63 


1.37 


0.6 


4.5 


0.70 


2.33 


HH 208B 


1.15 


3.90 


11.1 


3.87 


0.6 


4.5 


2.94 


13.0 


HH 208D 


0.85 


2.55 


9.96 


1.96 


0.6 


4.5 


2.31 


14.0 


HH 208E 


0.72 


2.28 


7.96 


1.84 


0.6 


4.5 


2.03 


10.8 


HH 208F 


0.25 


1.88 


7.96 


1.72 


0.6 


4.5 


1.90 


11.0 


HH 208J 


0.54 


2.28 


9.49 


1.84 


0.6 


4.5 


2.24 


13.5 


HM zUoiN 


3.31 


4.10 


9.53 


3.09 


0.6 


4.5 


1.50 


11. 1 


HH 208P 


0.02 


1.68 


7.66 


1.01 


0.6 


4.5 


2.38 


11.7 


HH 208R 


1.85 


2.55 


9.15 


1.60 


0.6 


4.5 


1.45 


13.0 


HH 208U 


0.0 


0.42 


1.82 


0.21 


0.6 


4.5 


0.63 


2.82 


128-248 


0.02 


1.28 


6.34 


1.01 


1.0 


8 


2.57 


13.6 


135-246 


0.33 


1.08 


3.20 


0.71 


0.6 


4.5 


1.10 


4.38 


137-239 


0.34 


1.48 


7.02 


0.95 


1.0 


8 


2.61 


15.4 


137-240 


0.33 


1.54 


7.32 


1.01 


1.0 


8 


2.76 


16.0 


140-239 


0.56 


2.48 


9.74 


2.37 


0.6 


4.5 


2.16 


13.0 


142-240 


0.54 


1.81 


5.28 


1.72 


0.6 


4.5 


1.43 


6.34 


143-225 


0.85 


3.90 


12.6 


5.15 


1.0 


8 


3.80 


32.4 


143-239 


0.51 


0.94 


1.49 


0.83 


0.6 


4.5 


0.54 


1.19 


144-237 


0.54 


1.20 


3.74 


0.89 


0.6 


4.5 


0.96 


4.98 


145-204 


0.18 


1.01 


4.17 


1.13 


1.0 


8 


1.31 


7.83 


152-229 


0.30 


2.15 


7.83 


3.15 


0.6 


4.5 


1.16 


8.51 


159-242 


0.69 


3.22 


12.0 


1.84 


0.6 


4.5 


3.88 


17.8 


161-246 


0.85 


2.48 


7.70 


1.54 


0.6 


4.5 


2.52 


10.8 



'^Extinction estimates from IChrvsostomou et al.l (|l997l ) and [Rosenthal et al.l (|200d ) - see text for more 
discussion. 



'Conversion to Ay based on lCardelli et al.l (jl989l ) 
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Fig. 2. — The estimation of H2 line brightnesses for 128-248. The boxes show the measured 
photometric points; errors are smaller than the boxes. The X's show the estimated continuum 
at the wavelengths of the H2 hues. The 2.12 /xm continuum was set equal to the 2.15 iim 
measurement. 
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The extinction corrections in this region are complex, and vary depending upon the 
distribution of intervening dust, which hes not only within the molecular cloud but also in 
the PDR and in foreground material. Extinction estimates are also complicated by reflec- 
tion off the back side of the nebula. Most of our objects, including two finge rs and several 



more compact structures, have some associated Fe II (1.64 /xm) emission (ISchultz et al. 



20071). For these objects, the 2.12 fim extinction can be estimated from the results of 



Chrysostomou et al.l (Il997l ). who used measurements o f the Fe II 1.257 and 1.644 /im transi- 



tions, and the extinction curve of ICardelli et al.l (119891 ). to obtain the extinction values of 0.6 
shown in column 6 of Table 1. For five regions, no Fe II emission is seen, possibly because 
of larger extinction. The extinction for these particular objects hasn't been estimated by 



Chrysostomou et al.l (Il997l ). or others. The best references are large beam e xtinction studies 



of the H2 Peak 1 source. The largest and most recent such study is that of iRosenthal et al. 



( 120001 ). based on the ISO measurement of 56 H2 transitions covering wavelengths of 2-17 ^m. 
They find Ak = 1.0 ± 0.1, a value we adopt here for 128-248, 137-239, 137-240, 143-225, 
and 145-204. For wavelengths other than 2.12 //m, the brightnesses have been extinction 
corrected using = A2.i2{ fii)"^'^^ (ICardelli et al.lll989l ). who also found that the shape of 
the extinction curve in the IR is independent of the value assumed for Ry. Usage of different 
plausible extinction corrections, within the uncertainties, does not appreciably alter our con- 
clusions. For reference, we show in column 7 the appr oximate visual extinc tions, Ay = 4.5 
and Av ~ 8, based upon Ry =Av /'E{b-v) = 3.1, and ICardelli et al.l (119891 ). The final 1.89 
/im H2 line-to-continuum ratios range from 0.25 to 1.1, with a median of 0.7. 



3. Shocked Emission Features: Line Ratio Analysis 



Early shock models of the Orion outflow invoke d planar C-type shock niqdels to explain 



the e mission from species such as H2 and CO (e.g. iDraine &: Robergdll982l . IChernoff et al. 



19821 ). C-type or "continuous" shocks occur at relatively low shock speeds {Vshock ^ 50 km s ^] 



in the presence of magnetic fields. A low ionization fraction allows ionized gas to cushion the 
shock in the neutral gas, limiting the neutral gas temperature to less than several thousand 
Kelvin and preventing significant dissociation. J-type or "jump" shocks generally occur 
at relatively high shock speeds {Vshock ^ 50kms~^) and usually dissociate molecular gas 
in the high temperature (T ~ 10'* — 10^ K) post-shock region. Molecules reform in the 
cooling post-shock gas at T~500 K. For pre-shock densities > lO^cm"^, H2 line ratios pro- 
duced in the reforming molecular gas may reach values higher than thermal values since 
H2 reforms in excited states , leading to a non-thermal cascade through rovibrational states 
jHoUenbach fc McKeelEgsgj ). 
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Observations of shocked H2O emission in Orion (e.g. lHarwit et al.lll998l ) seem to confirm 
tfie general picture that C-type shocks are responsible for the molecular emission in the 
outflow. Thes e studies converge o n preshock conditions n{B.2) ~ 10^ cm~^ and Khocfc ~ 
35kms"^ (e.g. IChernoff et al.lll982l ). It should be noted, however, that these studies fit data 
collected from a beam area covering an entire outflo w lobe 1' ~ 0-1 p c at Orion). Images of 



the shocked emission on sub-arcsecond scales (AB, IStolovy et al.lll998l ) show many emission 
features, each of which presumably has its own shock conditions. What is unique about the 
NICMOS observations of the inner region is that the 0."2 (~ 1.4 x 10^^ cm) resolution allows 
the isolation of individual shock fronts on the length scales expected for such shocks. These 
scales are expected to be ~ 10^^ cn a, depending on the presho ck density but with only a weak 
dependence on the shock velocity (IKaufman fc Neufeldlll996l ). This means that our deduced 
shock parameters are more likely to represent the local physical conditions, rather than an 
average over a number of shocks. 

H2 adaptive optics observations of the ambient mole cular cloud to the south-east of 
BN/K L have been carried out with 0."15 angular resolution ( Vannier et al.|[2001 : Kristensen et al, 
2OO3I ). The observed H2 1-0 S(l) brightness is well matched by C-shock models with shock 



velocities of 30kms ^ and pre-shock densities of 10 ^ cm 



but t he same models fall short 



of matching the 2-1 S(l) brightness by a factor ~2 fjVannier et al.ll200ll ). Higher shock ve- 
locity models improve the 2-1 S(l) brightness prediction, but provide a worse fit to the 1-0 
S(l) brightness. J-shock models produce the observed 2-1 S( l)/l-0 S(l) brightness ratio 
but h ave trouble reproducing the individual line brightnesses. iPineau des Forets fc Flower 
( 120011 ) suggested that non-stationary C-shocks can reproduce the high brightness and the 
large observed 2-1 S(l)/l-0 S(l) brightness ratios. However, non-sta tionary C-shock s have 
difficulty accounting for the proper motion velocity measurements of iDoi et al.l (120021 ). who 
found that over the velocity range of 20 - 400kms~^, different emission lines from the same 
object have similar velo cities. From comparison with the higher excitation H2 2-1 S(l) line, 
Kristensen et al.l ( l2003l ) found H2 clumps to have abrupt, south-facing edges exhibiting high 
excitation temperatures. Even C-shocks propagating into high densit y material can't account 
for the excitation temperature maxima and the line fluxes. This led Kristensen et al.l (|2003l ) 
to suggest that J-shocks propagating into material with a pre-shock density > lO^cm"^, 
plus an additional contribution from photo-dissociated mate rial, are required to explain 
these measurements. Emission from bow-shaped shock fronts (ISmith et al.lll99ll ). in which 
both C-shocks and J-shocks are responsible for the overall emission, may also help explain 
these observed brightnesses and line ratios. 

For our study, we use the ratio of the H2 line fluxes at 2.12 and 1.89 fim as a diagnostic 
of conditions in the emitting gas. The Camera 3 F190N filter transmission plot was exam- 
ined to determine the relative throughput of the 1-0 S(4) (1.892 /xm) and 2-1 S(6) (1.8947 
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/iin) rovibrational transitions. From the transmission plot, we estimate the NICMOS filter 
transmissions at the H2 wavelengths to be 64% and 98% respectively. Hence, we compare 
the measurements with model predictions for 0.64xB[l-0 S(4)]+0.98xB[2-l S(6)]. All three 
transitions require gas temperatures in excess of 1000 K in order to produce significant emis- 
sion. Such high temperatures imply that shock excitation is responsible for the H2 emission. 
A difficulty with interpreting the 2.12 and 1.89 /xm lines is that the transitions are relatively 
close together in excitation: the upper state of the S(l) line is at 6956K, and that of the 
S(4) and S(6) lines are at 9286K and 16880K, which can be a rather small baseline to fit. 

In an effort to characterize the shocked H2 emission, we have compared our observed line 
brightnesses and line ratios with standard models of shock emission: the C-shock model of 
Kaufman fc Neufeldl Jl996h and the J-shock model of iHoUenbach fc McKeel Jl989h . We first 
consider 128-248, a well-isolated source for which we obtain the highest signal-to-noise ratio 
at 1.89 and 2.12 jim. Because the large beam studies suggest shocked gas with densities near 
lO^cm-3, we computed the brightnesses of the 1-0 S(l) and 0.64xB[l-0 S(4)]+0.98xB[2-l 
S(6)] transitions in C-shocks with preshock densities of 10^, 10^ and lO^cm"^ and shock 
speeds up to 50kms~^; and in J-shocks with preshock densities of 10^, 10^ and lO^cm"^ 
and shock speeds up to lOOkms"^. In line with current practice, the models assume a) a 
planar shock, b) a magnetic field perpendicular to the direction of propagation of the shock, 
and c) the magnetic field str ength in microgauss equals the square root of the density in 



cm ^ (ITroland fc Heileslll986l ). The precise structure of the C-shocks depends on the field 



strength, as well as the ionization fraction, the grain size distribution, and the details of 
gas coolin g. We have previously exp lored the effects of varying these parameters on shock 
structure (IKaufman fc Neufeldlll996l ). We find that while the precise value of shock velocity 
and density determined from a line ratio may vary from those presented here, the range of 
intensities and line ratios possible in C-shock models is essentially limited by the temperature 
at which H2 dissociates. Thus our conclusion that C-shocks can explain the emission in most 
of the observed features is robust even if the precise shock parameters are different from those 
we have assumed. 

At this assumed magnetic field strength, low velocity shocks are C-shoc ks and not the 
lower magnetic field strength, non-dissociative, molecular J-shocks discussed by lWilgenbus et al 
(I2OOOI ). Typically, molecular J-shocks prod uce a factor of 10-1000 times fainter H2 emission 
in the lines discussed here than C-shocks ( IWilgenbus et al.ll2000l ). Since these line fluxes 
would be undetectable, we have not included non- dissociative, molecular J-shocks in our 
grid of models. Higher velocity shocks would be dissociative J-sh ocks, which we do consider. 
C-shocks are also known to be unstable to the Wardle instability (jWardle Ill990l ) arising from 
perturbations in the magnetic field direction, an effect which is not taken into account in the 
steady-state C-shock models presented here. The effects of this instability on the strengths 
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of H2 emission lines has been explored by iMac Low fc Smith I (119971 ) and iNeufeld fc Stone 



(119971 ). Both studies reached the conclusion that, for shocks over the range of densities 
and shock velocities we consider (i.e. shocks for which H2 is the dominant coolant), the 
instability has little effect on the predicted intensities of H2 lines. Thus our steady-state 
models should be sufficient for modeling the emission presented here. Clearly, changing the 
assumptions in the shock models - bowshocks instead of plane-parallel shocks, stronger or 
weaker magnetic fields, a different orientation of the magnetic field with respect to the shock, 
etc. - will result in different predictions for the H2 line strengths, perhaps resulting in better 
(or worse) agreement with our data. A complete exploration of this parameter space is be- 
yond the scope of this paper, but we can show that for some commonly-assumed magnetic 
field properties, plane-parallel C-shock models provide a reasonable fit to the data, while 
plane-parallel J-shock models generally do not. 

The predicted brightness ratio {0.64xB[l-0 S(4)]+0.98xB[2-l S(6)]}/B[l-0 S(l)] from 
each model calculation is shown in Figure 3 as a function of shock velocity. Also shown is 
the measured ratio, 0.19, for 128-248. The line ratio is consistent with either C- or J-shocks. 
However, the line ratio and brightness from this feature is best fit by a C-shock model with 
velocity 36kms~^ and preshock density of 6 x lO^cm"'^. The brightnesses, line ratio, and 
model predictions for 128-248 are listed in Table 2. The model values are well within the 
systematic measurement uncertainties. J-shocks produce absolute intensities which are too 
low to match the observed values. 
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Table 2. Observed Brightnesses and Model Fit for 128-248 



0.64xB[l-0 S(4)]+0.98xB[2-l S(6)] 


B[l-0 S(l)] 


0.64xB[l-0 S(4)]+0.98xB[2-l S(6)] 


[ergs~^ cm~^ sr"-*-] 




/B[l-0 S(l)] 


Observed 2.57 x lO'^ 


1.36 X 10-2 


0.190 


ModeP 2.60 x 10"^ 


1.31 X 10-2 


0.198 


^Model parameters: vs = 36kms~^ and n — 


6.3 X 10^ cm" 


which corresponds to an assumed 



magnetic field of 0.25 milli-Gauss. 
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Fig. 3.— Predicted {0.64xB[l-0 S(4)]+0.98xB[2-l S(6)]}/l-0 S(l) line ratio from C-shock 
and J-shock models with magnetic fields oriented perpendicular to the shock propagation 
direction and strengths in microgauss equal to the square root of the density in cm~^ 
( iTroland fc HeilesI Il986l ). Results are shown for C-shock models with n{B.2) = 10^, 10^ 
and lO^cm"^ and Vshock = 15 — 50kms~^, and J-shock models with 77.(H2) = 10^, 10^ and 
lO^cm"^ and Vghock = 30 — 100 km s^^. Also shown is the measured value of the line ratio 
for the Orion feature 128-248 (horizontal solid line), with the statistical uncertainty in the 
ratio indicated by the dashed lines. 
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The extinction-corrected line flux ratio for each of ten locations in HH 208 and thirteen 
other features is plotted versus their 2.12 fim H2 line brightness in Figure 4. The shock model 
curves are overlain. Within the systematic measurement uncertainties, all but one of the 
observed line ratios are consistent with C-shocks having preshock densities 10^ — 10^ cm~^ 
and shock velocities of 20 to 45 km s~^ . The narrow range of shock velocities is not surprising. 
Slower C-shocks produce much weaker H2 emission and would not have been detected. Faster 
C-shocks break down into J-shocks, again with much fainter H2 emission because the H2 is 
dissociated. The consistency with larger beam studies does suggest that these studies yield 
reasonable average shock parameters. Table 3 lists shock velocities and pre-shock densities 
for all the objects. These values were estimated by interpolating within the grid of calculated 
C-shock models at densities of 10^, 10^, 10^ and shock velocities of 20, 25, 30, 35, and 40 
kms~^, whose curves are plotted in Figure 4. Extrapolations for those objects just outside 
the grid were made using a few additional models with shock velocities of 45 kms~^. Two 
objects - 143-239 and HH208A - are clearly outside the C-shock grid. 

Examination of Figure 4 suggests that higher shock velocities may be correlated with 
lower pre-shock densities. A variety of statistical tests show that this correlation is significant 
at the 3-4(T level. The log of the pre-shock density falls by ~ 0.6 for each shock velocity 
increase of 10 kms~^. In order to search for structure such as the "hot edges" found by 



Kristensen et al.l ( 120031 ). we have constructed and examined images of the line flux ratio for 
all twenty three of the features included in Figure 4. We find no significant structure in the 
ratio images except for 145-204 and 159-242, where there is a peak in the line ratio offset 
by 1-2" from the peak H2 1-0 S(l) emission. 
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Fig. 4.— The log of the {0.64xB[l-0 S(4)]+0.98xB[2-l S(6)]}/[l-0 S(l)] brightness ratio 
plotted versus the log of the 1-0 S(l) brightness. Measured values for each of the shocked 
fingers with both 1.89 /xm and 2.12 fim emission are given by the pentagons. The symbols 
for the HH 208 positions are hollow. The error bars in the data are smaller than the 
data symbols. A trend arrow indicating the direction and magnitude of the change in 
the observed ratio and brightness if the 2.12 fim extinction was increased from 0.6 to 2.4 
mag, is given by the solid arrow. The symbol outlined by a box represents source 128-248. 
Model curves are shown for C-shocks (solid curves, labeled with log[n(H2)/cm~^] =4, 5, 6 
and shock velocity of 25-40 km s~^) and J-shocks (solid, dotted, and dashed curves, labeled 
with log[density/cm~'^] =5, 6, 7 and shock velocity of 30-100 km s~^). For both C- and J- 
shocks, the magnetic fields are oriented perpendicular to the shock propagat ion direction with 
strengths in microgauss equal to the square root of the density in cm~^ (iTroland fc Heiles 



19861 ). Factors contributing to uncertainties in the model predictions are discussed in §3. 
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Table 3. Estimated C-Shock Velocities and Pre-Shock Densities 



Position vs (kms ^)'^ log[nH2(cni ^)] 



TTTT Or\OT~) 

MM ZOoD 


42 


4.5 


TTTT Or\OT~\ 

MM 2UoU 


32 


4.8 


tin ZUofj 


36 


4.5 


TTTT OAOT7 

HH ZUob 


66 


4.D 


HH 208J 


32 


4.8 


HH 208N 


27 


5.2 


HH 208P 


40 


4.5 


HH 208R 


26 


5.5 


HH 208U 


42 


4.0 


128-248 


36 


4.8 


135-246 


45 


4.0 


137-239 


32 


4.8 


137-240 


33 


4.8 


140-239 


32 


4.7 


142-240 


42 


4.3 


143-225 


27 


5.8 


144-237 


38 


4.3 


145-204 


33 


4.6 


152-229 


27 


5.0 


159-242 


41 


4.5 


161-246 


42 


4.4 



^Assumes planar shock propagating into 
a perpendicular magnetic field with strength 
in microgauss eq ual to the square root of the 
density in cm~^ ( Troland &: Heileall986l ). 
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4. Results and Discussion 

In this section, we outline the morphological and emission characteristics of the H2 
features. We categorize the features into fingers - structures which either exhibit a clear bow 
shock morphology or features we believe to be bow shocks approaching us at a low angle - 
and knots - mostly including a variety of clumps in HH 208. 



4.1. Fingers 



The array of i n ner fi ngers which comprises the butterfly-shaped H2 emission flrst found 
by lBeckwith et al.l (Il978l ). extends over a 90" broad region. Most of the northern H2 flngers 
(AB) are outside of our field to the north, but there ar e additional fingers to the south, 
east of the Trapeziu i n (IMcCaughrean fc Mac Lowl 119971 ). The velocity measurements of 
Chrysostomou et al.l (119971 ) found that in addition to strong, broad H2 emission over the 
entire source, there are high velocity components confined to discrete condensations. The 
high velocity components are ascribed to additional 'bullets' similar to those imaged in the 
northern fingers by AB. At our higher angular resolution, the morphology of the inner fingers 
(Finger 1) is quite varied. Some of the objects are revealed to be bright, well-defined bow 
shocks; others do not display distinct bows of any kind. Most of the objects are found to 
have complex structure, with what appear to be internal shocks. 



4.1.1. 128-248 and 135-246 



128-248 is a bright, well-defined bow-shock in the southwestern portion of the finger 
array. Because it is so bright, and well-separate d from the main array by d ark dust lanes, 128- 
248 is an excelle nt candidate for many studi es. ISchultz &: Burton I (120071 ) found a FWHM of 
30 kms~ ^. while IChrvsostornou et al.l (119971. their object 1 0) fou nd it to have a FWZI of 150 



kms ^ . ISchultz fc Burton I (120071 ) and iGustafsson et al.l (l2003l . their region 11) found that 



its line profile has no secondary line peaks, which is suprising for a bow shock, and almost 
unique among the inner fingers. From its peak velocity of km s~^. IGustafsson et al.l (120031 ) 



concluded that 128-248 is moving in the plane of the sky, but were unable to determine 
in what direction. The H2 emission morphology is shown in Figure 5 - a classic bow-shock 
shape which appears to be moving to the southwest. From the line fiuxes, we deduce a shock 
velocity and pre-shock density of 36 kms~^ and 10^'^ cm~^ for 128-248. 



Fig. 6. — 2.12 /iin continuum plus H2 1-0 S(l) line map showing the locations of 128-248, 
135-246, 137-239, and 137-240 relative to the dust lanes and BN. The Trapezium is off the 
bottom of the image. North is up and east is to the left. 
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135-246 is a bright bow shock at the end of a faint finger emerging fro m a extended 



conti nuum "tail" of IRc4 (cf. Figure 1). It is prominent in H2 1-0 S(l), Fe II ( jSchultz et al. 



20071 ). and the F190N and F2 15N continuum baii dpasses. The object also appears in O I 



and possibly in S II emission (jSchultz et al.l 120071 ). Our H2 image is shown in Figure 5 - 
from the line ratios we deduce a shock velocity and pre-shock density of 45 kms~^ and 10^ *^ 
cm~^. The northeast boundary of the bright tip of the 135-246 bow shock is very sharp. 
We suggest this is probably due to the emergence of the finger from behind a region of high 
extinction: Figure 6 shows that the southwestern extent of the H2 finger is partially obscured 
by the same dark, curving dust lanes which are near 128-248 and the region south of BN. 



4.1.2. 137-239, 137-240, and 140-239 



The H2 emission from 137-239, 137-240, and 140-239 is shown in Figure 7. From our 
analysis, 137-239, 137-240, and 140-239 are all fit by shock velocities of 32 to 33 kms~^ and 
pre-sho ck densities of lO^-^" ^-^ cm'^. 140-239 is a small knot 3" south of IRc4. NICMOS H2 
images (jstolovy et al. 1998) suggest that this object may be a bo w shock, perhaps c oming 
towards us at a low angle because of its relative l y larg e blue-shift ( jSchultz et al.l 120071 ). The 
lower spatial resolution of IChrysostomou et al.l (119971 ) combined 137-239 and 137-240 into 
a single feature, which they identified as a high-velocity "bullet" (#5 in their list) with a 
FW ZI of 100 kms~^ . The two features were also observed together in H2 1-0 S(0) as object 
b of iLacombe et al.l (120041). The morphology i n that line was essentially identical to the 
morphology in 1-0 S(l). iGustafsson et al.l (120031 ) found these objects to have a peak velocity 
of -25 kms~^, althou gh they do not indi c ate w hether they were able to distinguish between 
the two bow-shocks. ISchultz fc Burton I (120071 ) find a velocity closer to -50 kms^^ for 137- 
239. From Figure 6, we suggest 137-239 and 137-240 may also be emerging from a dust 
lane, like 135-246. 
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Fig. 7.— H2 1-0 S(l) image of 137-239, 137-240, and 140-239. The maximum brightness 
in these features is 1.1x10"^ ergs s~^cm~^sr~^. North is up and east is to the left. 
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4.1.3. 142-240, 143-239, and 144-237 



These three objects are shown in H2 1-0 S(l) emission in Figure 8. 142-240 is a blunt 
bow-shaped object southeast of the star at the head o f IRc4. It is bright in H2 and Fe II, 



and much fainter in S II and O I (ISchultz et al.l 120071 ) . The Fe II emission ( ISchultz et al. 



20071 ) is more extended than the H2, suggesting these transitions sample different regions. 
For 142-240, we deduce a C shock velocity of 42 kms~^ and a pre-shock density of 10^'^ 
cm~^. 142-240 is accompanied on its eastern side by a fai nter, larger bow-sha ped region of 



Fe II emission - 143-239 - whic h is not seen in S II and O I (ISchultz et al.ll2007l ). From its H2 
emission, which is blue-shifted ( ISchultz fc Burton 1120071 ). 143-239 appears to arise from a J- 
shock: from Figure 4, the velocity and density would be in excess of 100 kms~^ and 10'' cm~'^ 
respectively. A deri sity this high might be expected to produce water masers, and indeed 
Gaume et al.l (119981 ) find a water maser within the er ror box of 143-239. 143-239 connects 
to 144-237 - a bright knot of Fe II (ISchultz et al.ll2007l ) and H2 emission 4. ^^4 to the northeast 



of 14 2-240. Similar to 143-239, 144-237 exhibits blue-shifted emission ( ISchultz fc Burton 



4.3 (^j^ 3_ 



20071 ). with a deduced shock velocity of 38 kms ^ and a pre-shock density of 10 
There is also faint Fe II emission slightly (~ 0."6) north of the H2 emission ( ISchultz et al. 



20071 ). If this Fe II is associated with 144-237, it may be that 144-237 is moving in the 
direction of BN or that the bow shock is asymmetric. 
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Rc4 

\ 


' 142-240 
144-237 143^239^ 



Fig. 8.— H2 1-0 S(l) image of 142-240, 143-239, and 144-237. The maximum brightness 
in these features is 5.5x10"^ ergs s^^cm^^sr"^. North is up and east is to the left. 
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4.1.4- 152-229 and 143-225 



Chrysostomou et al.l (119971 ) identified 152-229 as a bow shock - their bullet #6. Our 



H2, Fe II, and 2.15 /xm continuum emission is shown i n Figure 9. lO'Del 



et all (jl997aD 



found that S II emission in this region is blue-shifted. iGustafsson et al.l (120031 ) found a 
peak velocity of -21 kms^^ and a displacement of 0."2 between the emission peak and the 
location of the maximum velocity. From the direction of this displacement, they deduced 
that the shock is p ropagating towards BN ( roughly be t ween 2 and 3 o'clock in Figure 9). In 
disagreement with IGustafsson et al.l (120031 ). iDoi et al.l (120021 ) found from the proper motion 
of the S II emission, that 152-229 is m oving slightly north o f east - away from BN/IRc2 - 



with a transverse velocity of 50 kms~^. IChrysostomou et al.l (119971 ) found 152-229 to have 
a FWZl of 110 k ms~^. T he estimated shock velocity and pre-shock density are 27 kms~^ 
and 10^ ° cm~'^. iPaper ll noted that the H2 knot has a pointed cap of Fe II emission on 



the southeast side of the object (the blue arc in Figure 9), away from the putative exciting 
source (BN/IRc2) of the outfiow. The positioning and morphology strongly suggests a bow 
shock in which strong shocks producing the Fe II emission form on the leading surface of the 
bow while weaker shoc ks producing^ H2 emi ssion form behind it. The cap is also seen in the 
S II and I images of 



emission ( O'Dell et al. 



O'Del 



et al. 



(Il997al ) and is possibly also visible in high-velocity S II 



1997bl ): this may be the unlabelled S II knot northeast of 147-234. 



Chrysostomou et al.l (119971 ) also identified 143-225 as a bow shock - their bullet #8, 
with a a FWZI of 140 kms~^. 143-225 is among the most blue-shifted features in the 



outfiow (jSchultz fc Burton 1120071 ) and may be a bow shock approaching us at a low angle. 
The H2 emission is shown in Figure 10. The deduced shock velocity and pre-shock density 
are 27 kms~^ and 10^'^ cm~^. The shape suggests that it is a bow shock pointed slightly 
north- northwest, which would mean the origin of the feature would be somewhere to the 
south-southeast - roughly opposite to the direction of BN. 
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Fig. 9. — False color image of 152-229. Red is H2 1-0 S(l), green is the 2.15 /xm continum, 
and blue is Fe II. North is up and east is to the left. 
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Fig. 10.— H2 1-0 S(l) image of 143-225, 152-229, 159-242 and 161-246. The maximum 
brightness in these features is 1.8x10^^ ergs s~^cm~^sr~^. North is up and east is to the 
left. 
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4.2. Knots 



4.2.1. HH 208 



HH 208, approximately 7" west of BN, was fir st discovered by lAxon fc Tayloii (119841 ) . 
The S II and O I HST images of lO'Dell et al.l (Il997bl') clearly show a number of small features. 
Based on S II and O III images, lO'Dell et al.l (Il997al ) identified three knots in HH 208. The 
detection of optical features suggests that the extinction to HH 208 is lower than to other 
H2 features, implying it lies more in the foreground. Figure 11 shows our H2, Fe II, and 
continuum images, along with our knot identifications. The HH 208 H2 emission takes the 
form of discrete clumps, whereas the adjacent H2 emission has a more fingerlike appearance. 
It is difficult to discern what process has created this collection of features. 

Knot A is the faintest of the 2.12 //m H2 knots, but is the original HH 208 - seen in both 
S II and O I. lO'Dell et al.l (Il997bl ) showed images of knot A in several filters, and noted that 
a line drawn through HH 208 and HH 208NW terminates near the proplyd 154-240. They 
suggested that 154-240 may be the source of HH 208, but the line drawn (which is symmetric 
through HH 208NW but not through the rest of the object) also falls near IRc2 and radio 
sources "I" and "n". iDoi et al.l (|2002l ) found no net proper motion of knot A. The structure 



of the bright core of knot A did change in a disorganized fashion between 1995 and 2000, 
which corresponds to motions over a range of about 50 kms~^. Based on the high- velocity, 
blue-shifted emission lines, they further suggested that HH 208 is moving almost directly 
at us, rather than being connected to 154-240. In our data, knot A is the most extreme 
position - being beyond our grid of C-shock models. It is certainly higher velocity than any 
other HH 208 location, but may either be a high density (> 10^ cm~^) J-shock, or a low 
density (< 10^ cm~'^) C-shock. The position of knot B, serving as a bridge between knot A 
and the rest of the object, is suggestive of a relationship between the forbidden-line and H2 
emission. At 42 kms"^ and 10"^-^ cm , knot B has the second lowest density after knot U. 
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Fig. 11. — False color image of HH 208 with features discussed in this paper labelled. Red 
is H2 1-0 S(l), green is the 2.15 /im continum, and blue is Fe II 1.64 fim. The gold star 
symbols denote stars. North is up and east is to the left. 
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The arrangement of knots D-E-P-N-J-F is suggestive of the 



"ring" 



that Schild et al. 



( 119971 ) pointed out in the northeastern edge of the Orion H2 emission. The HH208 knots 
have deduced shock parameters in the 27-40 kms~^ and 10^'^ - 10^'^ cm~^ ranges, with no 
clear pattern. However, there may be some excitation gradient with distance from B. Around 
knots B-D-E, the purple features show where the H2 and the Fe II emission coincide and 
the purple disappears beyond E. Since forbidden line emission arises from fast J-shocks, in 
these regions at least, the spatial coincidence of H2 and Fe II emission may be inconsistent 
with our general result that the H2 emission in HH 2 08 arises in C-shock s. The knot of 
high- velocity S II emission designated HH 208NNW by lO'Dell et al.l (Il997al ) corresponds to 
the Fe II emission we find accompanying H2 kn ots D-E-F; th i s emiss ion can also be seen in 
combined S II and O I emission in Fi gure 2 of lO'Dell et al. (Il997bf) . Farther away, knots 



J and N together form Object #11 of IChrysostomou et al.l (119971 ). one of the regions from 
which they detected discrete high velocity H2 emission. Knot P shows neither forbidden line 
emission nor high-velocity H2 emission. 



Knots R and U together form HH 208NW (lO'Dell et al.lll997af ). Knot R shows opti- 



cal for bidden line emission, including high-velocity, blue-shifted S II emission ( lO'Dell et al. 
1997ah . The motio n of R (129-21 6) in the plane of the sky is 49 kms-^ and of U (126 



214) is 65 kms~^ (iDoi et al.l |2002| ). both roughly to the west on a path that would have 
recently traversed the B-D-E-N-J-F ring. The se proper motio ns put both these objects in 
the vicinity of BN/I/n around 1000 years ago (IDoi et al.ll2002l ). consistent with many other 
H2 features in the BN region, but significan tly earlier than the 500 year old BN/I/n break-up 
( iRodriguez et al.ll2005l : iGomez et al.ll2005l ). In our data, knot R shows the highest pre-shock 
density (10^'^ cm~^) but the lowest shock velocity (26 kms~^) in HH 208 - suggesting the 
H2 is not in the same region which produces the high velocity S II emission. Knot U ties 
with B for the second highest velocity (42 kms~^) and has the lowest density (10^'° cm~^). 
Although they may be unrelated to the B-D-E-N-J-F ring, knot U (or perhaps knot R) may 
be faster-moving material - a wind or knot - that impacted the ambient material to create 
the ring. Knot A could be a faster moving section of the expanding ring which is moving 
towards us. The ambient material would then have been a local H2 clump or even the core 
of a single, low mass star forming region. 



4.2.2. 159-242 and I6I-246 



159-242 and 161- 2 46 - in the southw e stern lobe of the outflow - are part of OMC Pk2 
dBeckwith et al.lll978h . ICustafsson et all J2003[ ) found 159-242 (their object 6) to have a 



peak velocity of +11 kms ^ and 161-246 (the western half of this knot is their object 19) to 
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have a peak ve locity of -15 kms ^. Applying shock models to their 2.12 fim 1-0 S(l) flux 



measurements, iVannier et al.l (120011 ) derived a pre-shock density of ~ 10^ cm~^, yielding 
a mass in 161-246 of 0.1 to 0.15 Mq, making it the most massive clump in their field. 
This led them to sugge st that the clump is a candidate site for low-mass star formation. 



Kristensen et al.l (120031 ) expanded upon that work by including 2-1 S(l) images; revising 



their shock models in light of t he new data led them to conclud e that the density is an 
order of magnitude greater than IVannier et al.l (120011 ) calculated. iKristensen et al.l (120031 ) 
also concluded that the 2-1 S(l)/l-0 S(l) flux ratio included a contribution due to radiative 
excitation from 6*^0 Ori, as well as to shocks. They were unable to reproduce both brightness 
and line ratios with a single type of shock, and therefore suggested that the shock contribution 
to the emission is composed of C-shocks in the interior of the clump, with J-shocks on the 
exterior. Our H2 emission is shown in Figure 10. From our analysis, 159-242 and 161-246 
are well fit by shock velocities of 41 and 42 kms~^ and much lower pre-shock densities of 10"^'^ 
and 10^'^ cm~^. This suggests that the material producing the 2.12 /xm 1-0 S(l) emission is 
insufficient to support even low mass star formation. We do find that the line ratio is about 
50% higher in a small region 0.9" east of the maximum H2 emission in 159-242, implying a 
hi gher shock velocity and lower density there, consistent with the exterior J-shock proposed 
by IKristensen et al.l (120031 ). 



5. Summary 

From 0."2 (90 AU) angular resolution HST NICMOS narrowband images of OMC-1, 
which resolve individual shocks, we estimate the brightnesses of H2 transitions at 1.89 and 
2.12 yum for 23 features. A comparison of the data with shock models shows that most 
of the data cannot be fitted by J-shocks, but are well matched by C-shocks with shock 
velocities in the range of 20-45 km s~^ and preshock densities of lO'* — 10^ cm~^. The narrow 
range of shock velocities is not surprising since both slower C-shocks and faster J-shocks 
produce weaker H2 emission and would not have been detected. Although there are many 
shock features in the OMC-1 region, most of the features appear to be well-characterized 
by a limited range of shock velocities and preshock densities, supporting the possibihty of 
a common origin. Additionally, these values confirm the findings of larger beam studies, 
which averaged over a number of individual shocks. Two objects - 143-239 and HH208A - 
are possibly due to J-shocks and the former does coincide with a known water maser. Optical 
forbidden line measurements of some features in HH 208 require fast J-shocks for excitation; 
we cannot explain this apparent discrepancy. 
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